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ABSTRACT 

Baryon interactions with bosonic dark matter are constrained by the potential for dark 
matter-rich neutron stars to collapse into black holes. We consider the effect of dark matter 
self-interactions and dark matter annihilation on these bounds, and treat the evolution of the 
black hole after formation. We show that, for non-annihilating dark matter, these bounds 
extend up to mx ~ 10^~^ GeV, depending on the strength of self-interactions. However, 
these bounds are completely unconstraining for annihilating bosonic dark matter with an 
annihilation cross-section of (cTav) ^ 10^^^ cm^/s. Dark matter decay does not significantly 
affect these bounds. We thus show that bosonic dark matter accessible to near-future direct 
detection experiments must participate in an annihilation or self-interaction process to avoid 
black hole collapse constraints from very old neutron stars. 
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1 Introduction 



Dark matter (DM) has been detected only via gravitational interactions. While there is 
overwhelming cosmological and astronomical evidence for new matter which could have a 
weak coupling to standard model fermions, the mass and couplings of this dark matter are 
not yet established. Asymmetric dark matter (ADM) has been proposed as a compelling 
framework to explain both the dark matter abundance and the baryon asymmetry [l -33 



It has been pointed out that models of asymmetric dark matter can be tightly constrained 



by the existence of old neutron stars 34 41 . The basic point is that dark matter will be 



captured by neutron stars due to DM-neutron scattering; if dark matter does not annihilate 
or decay, then it will continue to accumulate until it collapses into a black hole. The observation 
of old neutron stars can thus bound the dark matter-neutron scattering cross section (cr„x) for 
models with no dark matter annihilation or decay. It has also been demonstrated that bounds 



on non-annihilating dark matter shift with the introduction of self-interactions terms 42 -46 
In this work, we will study a more general question. We will focus on the range of anx, 
dark matter annihilation cross section {cav), dark matter decay rate (F), and self-interaction 
strength which could be consistent with observations of old neutron stars. 

As an initial point, we note that asymmetric dark matter is not necessarily non-annihilating. 
Asymmetric dark matter requires that the dark matter particle can be distinguished from the 
anti-particle, which implies the existence of a continuous unbroken symmetry under which 
the dark matter is charged. It is this symmetry which forbids the fermion Majorana mass 
term. However, asymmetric dark matter need not be the lightest particle charged under this 
continuous symmetry. For dark matter to be stable, it must be the lightest particle charged 
under some symmetry, but this may be a distinct Z2 symmetry. In this case, self-annihilation 
of the asymmetric dark matter would not be forbidden. 

Here we summarize the important features of the following analysis: 

I. Dark matter accumulation. The rate at which a neutron star captures dark matter (Cx) 
depends on the local dark matter density (px), the dark matter- neutron scattering cross 
section {(Jnx), and the dark matter mass {rrix)- The number of dark matter particles in 
a neutron star can be depleted by dark matter decay or annihilation. 

II. Black hole formation. Bosonic dark matter collected in a neutron star will form a black 
hole if the total energy is minimized at a radius which is less than the Schwarzschild 
radius. This condition is usually fulfilled by requiring that the dark matter be self- 
gravitating and exceed the Chandrasekhar limit. The Chandrasekhar limit grows with 
the strength of repulsive self-interactions. 

a. If the dark matter cannot thermalize quickly enough, it will not form a black hole 
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within the hfetime of the neutron star (see Appendix [B|). 

b. Thermahzed dark matter will collect within a radius rth, which determines the 
number of dark matter particles needed for the dark matter to self-gravitate. But 
if dark matter forms a Bose-Einstein condensate (BEC), then it will collect within a 
much smaller radius (rc) and fewer particles will be needed before the dark matter 
becomes self-gravitating and gravitational collapse occurs. 

Ill Destruction of the neutron star. The neutron star will be destroyed if the black hole 
grows large enough to consume the neutron star. The black hole will accrete baryonic 
and dark matter, but will emit Hawking radiation. To destroy the neutron star, the 
black hole must be large when it forms. 

a. If dark matter forms a BEC, then dark matter captured by the neutron star after the 
formation of a black hole will be efficiently accreted by the black hole, potentially 
compensating for the effect of Hawking radiation. 

b. If dark matter has a repulsive self-interaction, then the number of particles needed 
to form a black hole will increase. The resulting black hole, when formed, may then 
be large enough to continue growing. 

In section 2, we describe the accumulation of dark matter in neutron stars, including effects 
from dark matter decay and annihilation. In section 3, we describe the formation of a black 
hole from dark matter in neutron stars, including the effects of self-interactions. In section 
4, we describe evolution of a black hole which has formed in a neutron star, including the 
effects of baryonic and dark matter accretion and Hawking radiation. In section 5 we find the 
constraints on the parameter space of bosonic dark matter models from observations of old 
neutron stars. We conclude with a discussion of our results in section 6. 



2 Dark Matter Accumulation in a Neutron Star 



The dark matter capture rate (Cx) of neutron stars is given in [37,47 
(assuming T 



For mx < 0(10^) eV 



10^ K), dark matter within an old neutron star will be depleted by evaporation, 
and the rate of dark matter accumulation will also be suppressed. The probability for a dark 



matter particle to scatter while passing through a neutron star is given by 48 



P = 1 — exp 



- / VnO-nxdl 



(1) 



where r/n is the neutron number density and the integral is taken over the path of the dark 
matter particle through the neutron star. For small 0"„x, we have P ~ rjnAlcTnx- But P — )■ 1 
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for J rin(JnX<tt ^ 1. Taking a neutron star to have radius R = 10.6 km and baryonic density 
7.8 X lO^s GeV/cm^ |34||37||49) , this saturation occurs for ct„x > (r^at ~ 2.1 x 10"^^ cm^. 



For mx > GeV the capture rate is 37 



Cx ~ 2.3 X 10^^ Gyr-i (—) ( ^-^^^) f{anxmmx,mx,Vesc,v), (2) 

\n^x J GeV/cni'^/ 



while for mx < GeV, due to the effects of Pauh blocking, the capture rate is 37 



f{anx)(3{mx,mN,Vesc,v), (3) 



where px is the ambient density of dark matter. The factor f{crnx) is given by / = (Tnx/<^sat 
for anx < CTsat, and / = 1 for a„x > CTsat- The factor 

. l-exp[-6fe/^^)W(/i-l)^)] 

will take the value /3 ~ 1 for dark matter masses mx ^ 10^ GeV and for typical neutron star 
parameters, where /i = mx/friN, Vesc — 1-8 x 10^ km/s is the escape velocity from the surface 



of the neutron star, and v ~ 220 km/s is the ambient dark matter average velocity 34,37,50 

If tns is the lifetime of a neutron star, then the number of particles accumulated over 
that lifetime, Nacdtns), is determined by the dark matter capture rate, decay rate, and the 
rate of dark matter annihilation. We will first consider the case where dark matter does not 
annihilate, but does decay at rate F = r^^. In this case, the number of accumulated dark 
matter particles can be written as 

The tightest model-independent constraints on dark matter decay come from an analysis of the 
stability of dark matter halos. Measurements of halo mass-concentration and galaxy-cluster 
mass compared with simulations of dark matter halo mass distributions disturbed by dark 
matter decay constrain any dark matter lifetime to r > 10 Gyr for all velocities of the dark 



matter decay products 53 . Constraints on dark matter decay that heats the CMB are tighter 



for decay products with v > 0.6c 54 . We see that for r ~ 10 Gyr and a neutron star lifetime 
of tns ~ 10 Gyr, the number of accumulated dark matter particles is only suppressed by an 
0(1) factor. Because the minimum allowed dark matter lifetime is on the order of the lifetime 
of a neutron star, dark matter decay does not significantly alter the amount of accumulated 
dark matter, and thus does not significantly alter constraints arising from observations of 
neutron stars. 



3 



Henceforth, we will assume that dark matter does not decay. The accumulated number of 
dark matter particles can then be approximated as 55 , 56 



dNac 

dt 



Cx- 



2 

acc 



Vt 



th 



{cTaV) 



Tanh 



th 



(6) 



where Vth = (4/3)7rrj^ is the volume within which the dark matter is thermalized, here assumed 
to be of constant density, and {(Tav) is the annihilation cross section. If the effect of self- 



interactions is small, the thermalization radius can be written as 35,37,41 



Tth = 240 cm 



T GeV\ 



1/2 



105 K mx J 



(7) 



It is useful to determine the range of physical parameters for which the argument of the 
hyperbolic tangent in eq. [6]is greater than unity, indicating that the collection and annihilation 
of dark matter in the neutron star has reached an equilibrium. For a more complete treatment, 
see Appendix \K\ 



3 Black Hole Formation 

In order to form a black hole, the dark matter collected in a neutron star must become dense 
enough that the energy of the dark matter is minimized as the radius of the dark matter 
distribution approaches zero. If Ndm bosonic dark matter particles of mass mx are confined 
to a sphere of radius r, then the energy of a boson is approximately given by 

1 Gui^Ndm , 2'KGphmxr^ 

-c/ ~ \ ^ • («) 

r r 3 

The first term is the relativistic kinetic energy, and the second and third terms are the gravita- 
tional potential energy due to DM-DM interactions and DM-baryon interactions, respectively. 
The requirement of "self-gravitation" ensures the second term of eq. (jsj) is larger than the 
third, so the second term will dominate as the dark matter collapses. The Chandrasekhar 
limit then corresponds to the requirement that the second term dominate the kinetic term. 

However, this effective Chandrasekhar limit depends on the total local potential of the 
dark matter and is modified if dark matter has self- interactions. A A|(/)|^ term is generally 
not forbidden by any symmetry of the theory (in the absence of higher dimension terms. 



4 



stability of the potential would require this interaction to be repulsive, A > 0). With this 



self-interaction 35 , the number of self-gravitating particles required to form a black hole will 
be |5l||52| 

Of course this expression reduces to the simpler limit N chand ~ m^^/m?^ for non-interacting 
bosons when A = 0. Note, however that if A/327r > m\/m1^, then Nchand ~ X^^'^m^Jm^. As 
the A 1 01^ interaction is not forbidden by any symmetry, there is no reason to expect A to be 
very small. Thus, unless mx is quite large, one would expect the Chandrasekhar limit to be 
dominated by the interaction term. In this case, the Chandrasekhar limit on the number of 
particles is suppressed from the fermion case {N^l^J^'^ ~ m^Jm\) by a factor A^/^. 

The dark matter particles will become self-gravitating when their density exceeds that of 
the baryons in the neutron star. If the dark matter is confined to a region of radius r, then 
the number of dark matter particles required to achieve self-gravitation is given by 

N.-,(r) . '^^p., (10) 

where ph ~ 7.8 x 10'^® GeV/cm^ is taken as the baryon density in a neutron star. 



If dark matter forms a Bose-Einstein condensate 35 ,37 , then a large fraction of the dark 
matter will be confined to a radius which is much smaller than the thermalization radius of 
the dark matter. Thermalized bosonic matter at the core of a neutron star will form a BEC 
if the number of thermalized particles exceeds 

3\ fmxTV'^ /47rri\ , „,« / T 



The effect of self-interactions on the formation of a BEC is still not completely understood. 
In principle, self-interactions can affect the critical temperature and the size of the BEC state. 
A complete study of these effects is beyond the scope of this work. Instead, we will assume 
that the critical temperature and the size of the BEC state are unchanged by self-interactions 
of the magnitude which we will consider. If, due to self-interactions, dark matter does not 
form a BEC, then the analysis of Appendix [C] would be relevant. 

If the number of dark matter particles in the BEC phase is small, and the gravitational 
potential energy is dominated by the baryonic contribution, then the size of the BEC, Tc, 
can be approximated by equating the magnitude of the non-relativistic kinetic energy and 



gravitational potential energy 35 , 37 , yielding: 

1/4 /n^-\T\ 1/2 
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We have assumed that the contribution is small. Considering the ground state to have 
size Tc, the contribution of the term to the energy of the BEC state scales approximately 
as oc XN"^ /r^m\. This contribution is negligible for A <C 10~^^(mx/ GeV)^, implying that the 
value of the critical temperature and the size of the BEC state are essentially unchanged. From 



iVif f ^) = lO^^ ( —V^' . (13) 



eq. (|10|) we find that the number of particles in the BEC phase required for self-gravitation is 

mx J 

Assuming T = 10^ K, one finds Nchand > Ns-g if mx > 4 x 10"^^ GeV. For all mx of interest 
BECs will become self-gravitating well before they reach the Chandrasekhar limit. 

We may write the number of dark matter particles needed for black hole formation in the 
BEC phase as N^^^^^ms- We find 

^Blf/oL(^X,A,T) = N^Hand + NBEC (M) 

Note that, if dark matter forms a BEC, the neutron star must collect Nbec particles which lie 
within Tth and cause the formation of the BEC, as well as an additional Nchand particles which 
fall into a BEC of size and which collapse to form a black hole. For sufficiently large m^, 
one finds Ns_g{rth) > Nbec, in which case the dark matter collected within the thermalization 
radius will self-gravitate before enough dark matter is collected to form a BEC. This leads to 
the possibility that dark matter within the thermalized region will collapse without forming 



a BEC 46 , 57 . However, as dark matter in the thermalized region collapses it will also lose 



energy, which can result a lower temperature and higher density. This may subsequently lead 



to the formation of a BEC 58 



4 Hawking Radiation and Neutron Star Destruction 

The formation of a black hole within an old neutron star is not necessarily in conflict with 
observation - the black hole must also absorb the neutron star in a time much shorter than the 
lifetime of the neutron star, without flrst evaporating through Hawking radiation. We build 



on the analysis of 35 , 37 , but also consider the effect of dark matter accretion and repulsive 
self-interactions on the growth of the black hole. 

The evolution of the black hole's mass is governed by the equation 

dM,, Anp,{GM,f,f ^ rdM,f,\ 1 ^^^^ 



dt t;3 V J 153607r(GM6/,)2 
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where Mhh is the mass of the black hole and Vg is the sound speed of the neutron star (we take 



Vs/c ~ 0.1 [37]). The first term on the right hand side of eq. (15) is the Bondi accretion rate 
for baryonic matter, the second term is the rate at which the black hole accretes dark matter, 
and the last term is the Hawking radiation rate. {dMhh/dt)jjj^ will depend not only on how 
quickly the neutron star captures dark matter, but also on how quickly the black hole absorbs 
new dark matter captured by the neutron star. The initial black hole mass, Mf,hi is given by 

mxNs-g {Ng~g > Nchand) aud mxNchand {Ns-g < Nchand)- 

4.1 Rate of Black Hole Growth and Destruction 

If a black hole within a neutron star begins to grow, then the baryonic accretion rate will 
increase as M^^^, while the Hawking radiation rate will decrease as M^^^. In this case, we can 
approximate the time it will take for the black hole to consume the neutron star by assuming 
that the baryonic accretion rate dominates, neglecting dark matter accretion and Hawking 
radiation. We then find 

dt vl 
dM ~ AirpbiGMf 

/ 1 1 \ 



^ tnscoUa,se - (^^^^^ ^^^^ + M„ J ^ pj^,^, ' ^^^^ 

where M^s ~ 3.3 x 10^^ GeV is the mass of a heavy neutron star 59 . We consider the rate of 
collapse for a black hole with initial mass given by 

M--' ^ ^ 9.5 X 103' (^) (l . . (17) 



This initial mass yields a neutron star collapse time of 

A m 



= 2.6 X 10= yea.. (^) (l + 



PI 



327r m}-^ 



-1/2 



For all relevant regions of parameter space, this time of collapse will be small compared to the 
lifetime of an old neutron star. 

Similarly, if a black hole begins to shrink, the baryonic accretion rate will quickly become 
small, while the Hawking radiation rate will grow rapidly. Considering only the Hawking 
radiation rate, we find 

^ = -153607r(GM)2, 
dM V ; > 

= 51207rG2M,\,. (19) 
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We then find 

ti^g'^^ = 51207rG2f9.2xl03^-^V = 13Gyrf— v. (20) 
^ V rnx J \mx J 

We thus see that the black hole will evaporate in a time much shorter than the lifetime of the 
neutron star if nix ^ 1 GeV, for dark matter which forms a BEG. 



4.2 Black Hole Accretion of Dark Matter in the BEC Phase 



In order for dark matter to efficiently fall into a black hole after being captured by the neutron 
star, the impact parameter of the black hole must be small compared to the radius within which 
the dark matter settles In other words, additional dark matter will accumulate in the 
black hole at the rate it enters the neutron star if the size of the region where the dark matter 
particles settle is small compared to the black hole's impact parameter, bmfaii = 4:GMbh/ Voo- 
Here Voo is a dark matter particle velocity on approach to the black hole. 

After a black hole is formed from a BEG, Nbec dark matter particles will remain within 
radius r^^, and any further dark matter particles that collect in the star will fall into the BEG 
state. It can be shown that the radius of the BEG is smaller than the black hole's impact 
parameter [37], specifically hmfaii ~ 4rc, so in the case that all incoming dark matter settles 
into a BEG state, it will be efficiently captured by a black hole. Thus for dark matter which 
forms a BEG, the black hole dark matter accretion rate will equal the neutron star dark matter 
capture rate. 



bh 



dt 

The black hole will continue to grow if 
Cxmx > 

Cx > 



Cxmx. 



(21) 



DM 



4:TCpb{GmxNchandf 



153607r(G'mxiV, 



chand ) 



Gyr 



2.4 X 10 



36 



^ rnx \ 
\GeVJ 



1 + 



327r m'^ 



-1.5 X 10 



42 



GeVy 
mx J 



327r mj^ 



(22) 



If the black hole begins to grow, it will quickly absorb the entire neutron star. On the other 
hand, for the range of masses in which it is possible for Hawking radiation to dominate, the 
black hole will evaporate quickly. 
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A possible exception to even this bound arises if one assumes that Hawking radiation 
preferentially heats bosonic dark matter via dark-sector radiation. We refer to a detailed 



discussion of this effect in 37 . However, we note that, since a growing black hole will in fact 
grow rapidly, the Hawking radiation rate will quickly become small, implying that there will 
be very little heating of the dark or baryonic matter due to Hawking radiation. 

Bosonic dark matter with a large enough repulsive self-interaction cross section will have a 
larger mass at Chandrasekhar collapse and can avoid forming a small black hole that evaporates 
too quickly to destroy the neutron star. The result is an interesting phenomenon: very small 
repulsive self-interactions tighten neutron star collapse constraints on bosonic dark matter, 
but larger repulsive self-interactions loosen the same constraints. 



5 Bosonic Dark Matter Bounds From Neutron Star Col- 
lapse 

In this section we determine the constraints on cr„x for bosonic dark matter arising from the 
existence of old neutron stars, including the effects of dark matter self-interactions, dark matter 
annihilation, and dark matter accretion onto black holes. The exclusion contour bounds the 
region 

Nacc{0-nX,rnx,{(TaV),Px,tns,T) > NbH forms{mx , X,T) , 

dMsH , N 

> 0. (23) 



dt 



Figure [T] displays the exclusion contour in the (m^, plane if the dark matter can 

form a BEC, assuming that old neutron stars have lifetime t^s = 10 Gyr, core temperature 
T = 10^ K, and ambient dark matter density px = 10^ GeV/cm^ (this is an estimate for 
the dark matter density at the center of globular clusters [38|[39]). The various contours 
are for different choices of the self-interaction parameter A and the annihilation cross section 
(cxaf ) (thermally-averaged at the temperature of the neutron star). We note that an order of 
magnitude increase in (aav) corresponds to an order of magnitude relaxation of the bound, 
except in regions where the bound is significantly affected by dark matter accretion. Moreover, 
as the capture rate depends on px only through the factor px x cr„x, a scaling of the ambient 
dark matter density simply rescales the bound on anx- For example, if the dark matter density 
at the center of globular clusters is as small as 0.3 GeV/cm^, then the bound on anx would 
be weakened by a factor of ~ 3000. Note that observations of old neutron stars can only 
provide constraints in the region cr„x < o'sat.- For the entire relevant range of masses, bounds 
from neutron stars become completely unconstraining if (aav) x / > 10~^^ cm'^/s (thermally 
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mx[GeV] 

Figure 1: Neutron star collapse bounds for annihilating, self-interacting bosonic dark mat- 
ter that forms a Bose-Einstein condensate at globular cluster density px ~ 10^ GeV/cm^. 
From left to right the red, green, blue, and purple contours denote regions for which the self 
interaction parameter A = {0, 10~^°, 10~^^, 10~^^}, respectively. Sohd, dotted, dashed, and 
dot-dashed contours denote annihilation cross sections (o"af ) = {0, 10~^°, 10"''^, 10~^^}cm^/s, 
respectively. 



averaged at the temperature of the neutron star), because the neutron star can never capture 
enough dark matter for black hole collapse to occur. 

As previously noted, if the self-interaction term is increased (contours farther right in 
Figure [T]), the bound on high mass dark matter improves. As the Chandrasekhar bound 
increases with the self-interaction coupling, there is less Hawking radiation at the formation 
of higher mass black holes, and the black hole growth condition (eq. ([23])) is met for higher 
masses and lower scattering cross sections. 

Figure fl] also shows the excluded region (eq. 



29)) within which dark matter captured by 
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a neutron star will not thermalize, 



\GeVJ V T 



10^K\ flO Gyr 



tth 



(24) 



In the plot we assume a thermalization time scale of tth ~ Gyr. For rrix < 28 GeV, we 
reproduce the bounds of 37 . This analysis shows that old neutron stars in the center of 
globular clusters with a dark matter density px = 10^ GeV/cm^ 38 ,39 provide a bound 
on non-annihilating bosonic dark matter competitive with planned terrestrial direct detection 

10^ 



experiments for dark matter masses up to mx ~ 10 GeV [60, 61 . We note that future 
detection of neutron stars in regions of dark matter density larger than 10^ GeV/cm^ will 
result in appropriately rescaled bounds. 



5.1 Dark Matter Annihilation From a BEC State 

Thus far we have modeled dark matter annihilation as arising from a uniform distribution 
within the thermalization radius. However, if dark matter forms a BEC, then a significant 
fraction of the dark matter will be in the ground state. The BEC will be much denser; if dark 
matter can annihilate from this state, bounds on dark matter arising from observations of 
old neutron stars will be significantly weakened, if not entirely removed. When the BEC first 
forms, Tc/rth ~ 10"^. Thus, even a very small cross-section for dark matter to annihilate from 
the BEC state can result in a depletion of dark matter large enough to prevent a black hole 
from forming. Moreover, annihilation of dark matter in the BEC can heat the dark matter, 
also potentially obstructing black hole formation. 

However, one cannot determine the cross section for annihilation from the BEC state 
from {gav) in a model-independent way, since {(Tav) is determined by thermally-averaging 
the cross section at the temperature of the neutron star, T ~ 10^ K. Dark matter in the 
BEC state is initially much less energetic than dark matter in the thermalized region; the 
cross section for annihilation of dark matter in the ground state is thus model-dependent. 
As dark matter continues to accumulate in the BEC state, the dark matter will eventually 
become self-gravitating. Once the dark matter is self-gravitating, the size of the BEC will 
decrease as more particles fall into the ground state, causing the density of the BEC state 
to increase and causing the dark matter particles to have larger kinetic energy. By the time 
the Chandrasekhar bound is crossed, the dark matter will be relativistic. A more complete 
discussion of dark matter annihilation in the BEC state is beyond the scope of this work. 
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6 Conclusions 



We have studied the constraints that old neutron stars place on bosonic dark matter, allowing 
for self-interactions, decay, and self-annihilation of the bosonic dark matter. Observations 
of old neutron stars imply bosonic dark matter with a mass ~ kev - 10^ GeV detected at 
terrestrial experiments will have a minimum annihilation or self-interaction term. For example, 
we show that a neutron star of age tns = 10 Gyr found in a globular cluster with dark matter 
density px = 10^ GeV/cm^ will not constrain bosonic dark matter if the dark matter has an 
annihilation cross section (aav) > 10~^^ cm^/s (thermally-averaged at the temperature of the 
star, T ~ 10^ K). If dark matter has even a small cross section to annihilate from a BEG state, 
then constraints from neutron stars can be weakened even more. These bounds are thus most 
relevant if the unbroken symmetry which stabilizes an asymmetric dark matter candidate is 
continuous; if it is broken to a Z2 symmetry (even weakly), then self-annihilation is permitted 
and these bounds can be weakened considerably. Gonversely we demonstrate that permitted 
dark matter decay, which is constrained by the evolution of dark halos, will not significantly 
relax neutron star bounds on bosonic dark matter. 

We also show that even small self-interaction terms can dramatically weaken bounds on 
asymmetric dark matter. These bounds are thus most constraining if there exists some (at least 
approximate) symmetry which can suppress a quartic self-interaction term. However, very 
small self-interactions can result in even more constraining bounds, by causing the formation 
of larger black holes which grow rapidly. 

It is interesting to note that these bounds can extend up to large mx- For such large masses, 
one cannot easily tie the dark matter asymmetry to the baryon asymmetry. Nevertheless, 
asymmetric dark matter with a small (or vanishing) annihilation cross section provides an 
interesting candidate for non-thermal dark matter, and observations of old neutron stars can 
provide significant constraints on these models. 

Note added: While this paper was being completed, fSS] appeared, which also discusses the 
effect of dark matter accretion on the evolution of a black hole. In 58 , it is also argued that 
neutron star observations cannot bound dark matter particles with a large mass. Note that this 
argument refers to potential bounds on dark matter which self-gravitates in the thermalized 
region before forming a BEG. The bounds we have described for BEG dark matter at high 
mass are not affected by this argument. These bounds instead arise from a consideration of 
the effect of self-interactions, and the effect on the black hole's evolution of the accretion of 
dark matter in the BEG phase. 
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A Dark Matter Capture and Annihilation Equilibrium 

It is useful to determine the range of parameters for which the total dark matter in the neutron 
star will reach an equilibrium, at which point dark matter will annihilate at the same rate it 
collects in the neutron star. This happens when the argument of the hyperbolic tangent in 
eq. ^ is greater than unity - in this range the hyperbolic function will evaluate to unity and 
the formula for the number of accumulated dark matter particles simplifies considerably. For 
mx > GeV, 



(3 



X 1 10-45 cni3/s ■ 103 GeV/cm3 ' ^^^"^^ j ' 



and for mx < GeV, 



10 Gyr 



/3 



1/2 



1^10-45 cm3/s ■ 103 GeV/cm3 ' ^^"""^V ' 

In figure [2| we plot the region in the (mx, (o"af) x /)-plane such that the neutron star is in 
equilibrium, assuming tns = 10 Gyr, T = 10^ K and px = 10^ GeV/cm3. Even for small 
values of ct„x and {(Tav) (including any values which could be detected by current and planned 
observations), we may use the approximation 



Nacc ~ VCxVth/{aaV) 

/47rCx(240 cm)3 / T GeVV^^ 



3{aav) yiO^K mx J 
But if dark matter has not reached equilibrium, we instead find 

Nacc ~ Cxtns- (2 



(27) 
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Figure 2: In the dark shaded region of the {mx, {cTav) x /)-plane, the neutron star is in 
equihbrium, assuming tns = 10 Gyr, T = 10^ K and px = GeV/cm^. 



B Thermalization Time 

For the dark matter particles to achieve self-gravitation, the dark matter must thermalize with 
the neutron star on a time scale comparable to the neutron star lifetime. For mx ^ GeV, the 



thermalization time is 37 



tth = 5.4 X 



10 years (—j j^— 




(29) 



where / = Onxjosat. if o"nX < o"sat., and / = 1 otherwise. Here, Onx is the cross section for 
DM- neutron interactions, Ogat. ~ 2.1 x 10~^^ cm^ and T is the core temperature of the neutron 



star. Thus eq. (29) presents a firm exception to the application of the neutron star collapse 



bound (eq. (23)): if the dark matter does not thermalize over the lifetime of the star, it will 
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not form a black hole. 

It is important that the saturation of the likelihood of dark matter scattering implies a 
maximum dark matter mass for which neutron star bounds are applicable; beyond this mass 
the dark matter will not thermalize: 

/ rp \ 1/2 / , \ 1/2 
(max) _ T . ^ TnTn^M / \ / 



= 1.4 X lO'GeV j (^^j . (30) 

C Dark Matter Which Does Not Form a BEC 

If dark matter does not form a BEC, then the analysis is slightly different. The number of 
dark matter particles in the thermalized region required to achieve self-gravitation is given by 

For most regions of interest, dark matter which does not form a BEC will collapse when 
it becomes self-gravitating. But if A is large enough, thermalized boson distributions will 
collapse only when they reach the Chandrasekhar limit. In figure [3} we plot the region where 
Ng^^g > Nchand in the {mx, A)-plane for T = 10^ K, in the case where a BEC does not form. 

The number of particles required for a black hole to form in the case where a BEC does 
not form is then given by: 

^BHforms{f^X,\T) = max[Nchand, N^-g]- (32) 

If Nchand > ^s-g, then the initial mass of the black hole is the same as the value given in 
eq. 17 But if Ng-g > Nchand, then the initial black hole mass is given by 

mIT""^ = mxNt'r^ = 4.8 X 10^^ GeV (^^) . (33) 

In this case, if the black hole grows, the time required for the neutron star to be destroyed is 

CX» = 5.1 X 10- years (^)''^(1^)''\ (34) 

and the black hole will quickly destroy the neutron star. If the black hole evaporates, the time 
required to complete the evaporation process is given by 

= 51207rG= Us X lO* l^'i^ 
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Figure 3: Parameter space in the (mx,A) plane where Ng^^g' > Nchand, if a BEC does not 
form. 



If dark matter does not form a BEC, then the black hole will evaporation time will be short 
if mx > 10^ GeV. 

The dark matter thermalization radius is much larger than the black hole impact param- 
eter; thus dark matter will not be efficiently captured by a black hole for dark matter that 
does not form a BEC. Instead, the black hole will continue to capture dark matter through 
Bondi accretion. The dark matter accretion rate can then be written in terms of the coupled 
differential equations. 



dM, 



bh 



dt 



SmxNDMiGMi 



bh 



DM 
dNpM 

dt 



Cx- 



^srth 
1 /dM, 



bh 



(36) 



mx \ dt y 

However in the case of non-BEC dark matter black holes, the Hawking radiation rate 
exceeds the baryonic Bondi accretion rate only for 

T 



mx > 5.8 X 10" GeV 



105 K 



(37) 
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To accrete enough dark matter for the dark matter Bondi accretion rate to be of the same 
order of magnitude as the baryonic Bondi accretion rate would require about as much time 
as was required for the black hole to form in the first place. But as we have seen, for the 
range of masses where the Hawking radiation rate dominates when the black hole is formed, 
the black hole will evaporate relatively quickly, before the dark matter accretion rate becomes 
appreciable. We thus find that, for the case where dark matter does not form a BEC, we may 
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Figure 4: Neutron star collapse bounds for annihilating, self-interacting bosonic dark mat- 
ter that does not form a Bose-Einstein condensate. From left to right the red, purple, and 
orange contours denote regions for which the self interaction parameter A = {0, 10~^^, 10~^}, 
respectively. Solid, dotted, dashed, and dot-dashed contours denote annihilation cross sections 
(a^z;) = {0, 10-5°, 10-^^ 10-^2}(,jn3/g^ respectively. 

essentially ignore the effect of dark matter accretion on the evolution of the black hole and 
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the condition for a black hole to destroy the neutron star becomes, 



dt 



m^.- (38) 



3 

MBH=MBHi "^S 



153607r(GmxMax[iV,;,,„,, Nr!,]f 

In Figure |4| we plot exclusion contours in the (m^ , <Jnx ) plane if the dark matter cannot 
form a BEC, assuming that old neutron stars have lifetime tns = 10 Gyr and core temperature 
T = 10^ K, and assuming an ambient dark matter density of px = 10^ GeV/cm^. As expected 
from Figure [3| the self-interaction coupling does not affect the bounds substantially until 
A > 10^^^ - at this value the number of particles required for self-gravitation is exceeded by 
the number of particles required by the Chandrasekhar bound. 
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